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Abstract 

We calculate the transverse and time-time components of the instantaneous gluon propagator in 
Coulomb gauge QCD by using an SU(3) quenched lattice simulation on isotropic and anisotropic 
lattices. We find that the gluon propagators suffer from strong discretization effects on the isotropic 
lattice; on the other hand, those on the anisotropic lattices give a better scaling. Moreover, on 
these two type of lattices the transverse parts are significantly suppressed in the infrared region 
and have a turnover at about 500 [MeV] . The high resolution to the temporal direction due to the 
anisotropy yields small discretization errors for the time-time gluon propagators, which also show 
an infrared enhancement as expected in the Gribov-Zwanziger confinement scenario. 
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I. INTRODUCTION 



The Coulomb gauge with no negative metric provides a very clear physical picture in the 
sense that the color-Gauss's law can be formally solved, only transverse degrees of freedom 
appear as dynamical degrees of freedom, and Fock space is well defined. The prominent 
feature taking the Coulomb gauge is that an instantaneous interaction, which is requisite 
for color confinement, shows up in the Hamiltonian. In the Gribov-Zwanziger scenario, the 
path integral is dominated by the configurations near the Gribov horizon where the lowest 
eigenvalue of the Faddeev-Popov (FP) ghost operator vanishes [1]. The lattice simulations 
show the enhancement of near-zero modes of the FP eigenvalues [2, 3]. Accordingly, the 
color-Coulomb instantaneous interaction bears a confining force. It has been also confirmed 
by the lattice simulations that the color-Coulomb potential rises linearly at large distances 
and its string tension is larger than the string tension of the static Wilson potential [4- 
7], which is expected from the Zwanziger's inequality [8]. In addition, the color-Coulomb 
potential can be reevaluated by inverting the FP ghost matrix and this analysis has shown 
that its string tension almost saturates the Wilson string tension [9] . 

Exploring the gluon propagator is a central issue for studying the confinement mechanism 
in QCD. The transverse would-be physical gluon propagator is expected to be suppressed 
in the infrared (IR) region due to the proximity of the Gribov region in the IR direction 
in the Gribov-Zwanziger scenario [10]. There have been a lot of lattice studies and func- 
tional analyses in the Landau gauge (see, for instance, Ref. [11, 12] and references therein). 
Contrastingly, there are few lattice studies on the instantaneous gluon propagator in the 
Coulomb gauge [13-16]. In the continuum theory, on the other hand, there are many works 
on the variational approach in the Coulomb gauge, which are very useful to study the spec- 
trum of the hadronic bound states and the Green's functions [17-22], in addition to the 
functional analysis [23-25]. 

The Coulomb gauge fixing condition is imposed on each time slice and it does not intro- 
duce correlations between neighboring time slices. Accordingly, the equal-time correlation 
between gauge fields at different points 

(A*(x,t)Al(y,t)) (1) 

is well defined in the Coulomb gauge. Recent lattice studies of the instantaneous gluon 
propagator revealed that it shows scaling violation [15, 16]; namely, the gluon propagator 
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calculated at different lattice couplings do not fall on top of one curve after multiplicative 
renor malizat ion . 

In order to circumvent the problem of scaling violation, the authors of Ref. [15] have ex- 
tracted the instantaneous gluon propagator by eliminating the dependence of the unequal- 
time propagator. It has been concluded that the instantaneous transverse gluon propagator 
D tT (p) is multiplicatively renormalizable in the Hamiltonian limit and the numerical data of 
it are well fitted with the Gribov-type form of the propagator 1 



where M is a fitting parameter at which the propagator shows a turnover. The method 
was applied to the transverse gluon propagator successfully both in 2+1 and 3+1 dimen- 
sional SU(2) Yang-Mills theory [15, 26] while it does not improve the scaling violation of 
the temporal gluon propagator since the time-time component of the unequal-time gluon 
propagator is energy independent even after the residual gauge fixing [27]. 

As another approach, a new momentum cut is introduced in Ref. [16] in addition to the 
cone cut and the cylinder cut. High momentum data that suffer from discretization errors are 
excluded from the analysis of the instantaneous propagators by this new cut. Combined with 
the matching analysis given in Ref. [28], it has been shown that this procedure successfully 
reduces the scaling violation for the transverse gluon propagator, but it fails for the time-time 
component of the gluon propagator. 

In this study, we take a different route to show that the instantaneous gluon propagator 
is multiplicatively renormalizable. The problem of scaling violation of the instantaneous 
propagator can be seen even at the tree level with a finite temporal lattice spacing as 
was discussed in Ref. [15] (we shall briefly review the point in Appendix B). The reason is 
that the energy integral does not run from — oo to oo but from —Tr/a T to ir/a T on a finite 
lattice, and this introduces the spurious \p\ dependence on the free instantaneous propagator. 
Therefore, we expect that the instantaneous propagator is multiplicatively renormalizable 
in the Hamiltonian limit a T — > 0. To make this point clear, we calculate the transverse and 
the time-time components of the instantaneous gluon propagator on anisotropic lattices and 

1 In this paper, we use the same symbol for the instantaneous propagator and the unequal-time propagator, 
but the reader may distinguish them by their argument: the instantaneous propagator does not depend 
on X4 (or p4 in momentum space) but the unequal-time one does. 
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show how the scaling violation becomes milder as the anisotropy increases, i.e., as we get 
close to the Hamiltonian limit. 

The organization of this paper is as follows. In the subsequent sections, we describe the 
lattice observables, the space and time components of the instantaneous gluon propagator, 
and the lattice setup of our numerical simulations. In Sees. IV and V, the numerical results 
of the transverse and the temporal components of the instantaneous propagators on the 
isotropic lattice are reported. Section VI is devoted to show the results for the transverse 
propagator on the anisotropic lattices. In the subsequent Sees. VII, we discuss the IR and 
the ultraviolet (UV) behavior of the propagator by making the power law fitting. The 
anisotropic lattice results for the temporal gluon propagator are given in Sec. VIII, and the 
IR and the UV fittings are examined in Sec. IX. We attempt to extract the color-Coulomb 
string tension from the instantaneous temporal gluon propagator in position space in Sec. X. 
The conclusions are drawn in Sec. XI. 



II. INSTANTANEOUS GLUON PROPAGATOR 

We calculate the transverse and the time-time components of the instantaneous gluon 
propagator, 

Df u (x-y) = (Al(x,t)A b u (y,t)), (3) 



in the momentum space, 

DfM = ^jfp (5) 
where the gauge fields are related to the link variables through 

4 at (*,t) 



^=^^,-f§)D tr (p), (4) 

cab ^44 (P) 



llaU ^ U,(x,t)-Ul(x,t) 



2ig a^ 



(6) 

traceless 



The instantaneous gluon propagator is evaluated on each time slice and we average over all 
time slices. The lattice momenta fc M are discretized and take integer values in the range 
(— Lft/2, L^/2]. The lattice momenta and the continuum ones are related via 

2 f-nk,, 



P„ = — sin — , (7) 



where a M (L At ) are a a {Lo) for // = 1 to 3 and a T (L T ) for \x = 4, respectively. 

We note that the unequal-time propagator D^ v has mass dimension 2 in momentum 
space, while the instantaneous one has mass dimension 1. This is because the instantaneous 
propagator is obtained by integrating the unequal-time propagator over the time component 
of the four momentum; 

D,AP) = I ^D, v (p,p,). (8) 

The unrenormalized transverse gluon propagator D\l t , which is measured by lattice simu- 
lations, is related to the renormalized propagator via the multiplicative renormalization, 

Dr(f, y) = a a Z tx (a a , v)D{l t (pa a ), (9) 

where /i is the renormalization point. We expect that the renormalized propagator is inde- 
pendent of the lattice spacing in the scaling regime. As we shall see later, the multiplicative 
renormalizability does not hold for finite temporal lattice spacing and we have to take the 
Hamiltonian limit, a T — > 0. 

The color-Coulomb potential plays a crucial role in the Coulomb gauge QCD. It was 
shown that the time-time component of the gluon propagator can be decomposed into the 
instantaneous part and the noninstantaneous part [29], 

D u (x,t) = V c (x)S(t)+P(x,t). (10) 

The first term in the right-hand side represents the color-Coulomb potential, which is defined 
as the vacuum expectation value of the kernel of the instantaneous interaction, 

V c (x-y)6 ab = ({M-\A]{-dt)M-\A])f s ). (11) 

Here M~ l is the Green's function of the Faddeev-Popov ghost operator. In Eq. (10), P(x, t) is 
assumed to be nonsingular at t = as opposed to the first term. It has been shown that both 
V c and P are renormalization-group invariant [30] ; namely, the renormalization constant for 
the color-Coulomb potential can be set to be 1 in the continuum limit. We expect that 
the color-Coulomb potential V c can be extracted from the instantaneous temporal gluon 
propagator as 

V c (x) = a T D M (x), (12) 

when we are in the scaling region and the lattice spacing is small enough, where a T comes 
from the 5 function. 
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III. LATTICE SETUP 



The lattice configurations are generated by the heat-bath Monte Carlo technique with 
the standard Wilson plaquette action, 

S= 7^ 2 ^Tr(l-^(n))+/3e B ^^eTr(l-C/ i4 (n)). (13) 

?T..i<j <3 n,i<3 

Here U^ u {n) indicates the plaquette operator, and (3 = 2N c /gQ is the lattice coupling. On 
the isotropic lattice, the bare anisotropy is 1 and the action can be written in a familiar 
form 

S = (3 ^eTr(l-fVH). (14) 

n,/i<i> 

The renormalized anisotropy £ is defined as the ratio of the spatial lattice spacing to the 
temporal lattice spacing, £ = a a /a T . The ratio of and £ can be determined nonpertur- 
batively by matching the spatial and the temporal Wilson loops on anisotropic lattices. We 
use the relation obtained by Klassen for the range 1 < < 6 and 5.5 < (3 < oo [31]. We 
adopt the values of the lattice spacing given in Ref. [32] for £ = 2 and in Ref. [33] for £ = 4, 
where the static quark potential was measured to set the scale. For the isotropic lattice, the 
scale is set by using the scaling relation obtained by Necco and Sommer, with the Sommer 
scale parameter r$ = 0.5 [fm], which is applicable in the range 5.7 < /3 < 6.92 [34]. 

In our simulations, the first 5000 sweeps are discarded for thermalization, and we mea- 
sured the instantaneous gluon propagator for 50 — 100 configurations, each of which is 
separated by 100 sweeps. All the lattice parameters are given in Table I. 

In the Coulomb gauge the transversality condition 

d i A l {x,t)=0 (15) 

is imposed on the gauge fields at each time slice, where i runs from 1 to 3. On a lattice, 
gauge configurations satisfying the Coulomb gauge condition can be obtained by minimizing 
the functional 

Fu ^\ = 7?Z> eTr ( X " l9\Z,t)Ui(£,t)g(£+T,t)) (16) 

with respect to the gauge transformation g(x,t) e SU(3) on each time slice. The functional 
derivative of Eq. (16) with respect to g yields the transversality condition VjA- at (x, t) = 0, 
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TABLE I: The lattice couplings, the spatial and the temporal lattice extents, the lattice spac- 
ings, the lattice volumes in physical units, and the number of configurations used to evaluate the 
instantaneous propagators. 



£ = a a /a T L a L T /3 a a x [GeV] a a [fm] V^fm 4 ] Number of configurations 
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where V is the lattice backward difference, and it reproduces the Coulomb gauge condition 
in the continuum limit. The Coulomb gauge fixing has been done using an iterative method 
with the Fourier acceleration [35], and we stop the iterative gauge fixing if the violation of 
the transversality becomes less than 10~ 14 ; 

* = {N 2 l _ l)L { £^^^t)f < 10- 14 . (17) 

This stopping criterion is applied for each time slice. We note that the accuracy of the gauge 
fixing is crucial for the transverse propagator to see the IR suppression, which is discussed 
in Appendix A. 

In order to reduce lattice artifacts, we apply the cone and cylinder cuts to the momenta 
[28]. The cone cut is necessary to address finite volume effects that are seen in small 
momentum data. On the other hand, the cylinder cut reduces artifacts due to the broken 
rotational symmetry on lattice. The statistical errors are estimated by the jackknife method. 



IV. INSTANTANEOUS TRANSVERSE GLUON PROPAGATOR ON THE 
ISOTROPIC LATTICE 

Although the problem of the scaling violation for the instantaneous gluon propagator has 
already been discussed in Refs. [15, 16], we here show the lattice result for the instantaneous 
transverse gluon propagator on the isotropic lattice to clarify the issues. The instantaneous 
transverse gluon propagator on the isotropic lattice at (3=5.7 and 6.0 is drawn in Fig. 1. 
The propagator is normalized such that D tr (|p| = 2[GeV]) = 1. 

We observe that D tr has a maximum at p = 0.4 — 0.5 [GeV] irrespective of the lattice 
coupling and it decreases with the momentum in the IR region. This is a striking feature 
of the transverse gluon propagator. The instantaneous propagator is defined as the energy 
integral of the unequal-time propagator, 

D(\P\) = J ^D(?,P*)- (18) 

For a massless particle and a massive particle, -D(|p|) = l/(2|p|) and D(\p\) = 
1/(2a/|p| 2 + m 2 ), respectively. Thus the instantaneous transverse propagator can be in- 
terpreted as the inverse of the energy dispersion relation of the would-be physical gluons. It 
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FIG. 1: The instantaneous transverse gluon propagator in physical units at (3 = 5.7 and /3 = 6.0. 
The propagator is renormalized at \p\=2 [GeV]. 

implies that the propagator at vanishing momentum corresponds to the inverse of the effec- 
tive mass of the gluon. The IR suppression of the instantaneous transverse gluon propagator 
means that the gluons have momentum dependent effective mass M(p) and it diverges in 
the IR limit, lim^o M(p) = oo, indicating the confinement of gluons. 

In addition to the bump structure of the transverse propagator, the two curves corre- 
sponding to different lattice couplings cross at the renormalization point and deviate from 
each other in the small and the large momentum regions. This is not seen in the Landau 
gauge gluon propagator (see e.g. Ref. [28]). Such a behavior of the propagator casts doubt 
the validity of the multiplicative renormalizability for the instantaneous gluon propagator. 

Taking a closer look at the raw results of the numerical simulations gives us a clue to 
cure scaling violation. Assuming multiplicative renormalization for the propagator, the 
renormalized dressing function of the transverse gluon propagator, = \p\D^, is 
related to the bare unrenormalized dressing function Z\l t via 

Z£(p;fi) = Z(a (T ^)Zf at (\p\a (T ) (19) 

in the scaling regime. Here Z(a a ,fi) is a renormalization constant. It is easily read off from 
this relation that in the log-log plot of the dressing function of the propagator, converting 
from lattice units to physical units corresponds to the parallel displacement in the horizontal 
direction and the renormalization of Z tr corresponds to that in the vertical direction. If 
the propagator is multiplicatively renormalizable, the different curves associated with the 
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different lattice couplings can fall on top of each other by a parallel shift of the curves in 
the horizontal and the vertical directions in a double-log plot. 

The left panel of Fig. 2 shows the dressing function of the unrenormalized transverse gluon 
propagator in lattice units. It is apparent that the two curves cannot coincide by adopting 
any scaling relation or by imposing any renormalization condition since such manipulations 
correspond to the horizontal and the vertical shifts in the log-log plot of the dressing function 
but not the rotation or deformation of the curves. Therefore, the scaling violation of the 
transverse gluon propagator is purely due to discretization errors. 

Furthermore, we observe that the IR behavior of the dressing function at different cou- 
plings shows the same behavior, while in the UV region the slope of the curves differs. The 
right panel of Fig. 2, in which the dressing function renormalized at \p\=l [GeV] is plotted 
as a function of the physical momentum, illuminates such a tendency; the two curves almost 
fall on top of each other in the IR region while the deviation between them is pronounced 
in the UV region. This indicates that the scaling problem of D tr resides in the lattice data 
in the UV region. 
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FIG. 2: The dressing function of the unrenormalized instantaneous transverse gluon propagator, 
Z\l t {p) = \p\Dll t (p), at different lattice couplings in lattice units (left). The dressing function 
renormalized at |fj|=l [GeV] in physical units (right). 

This kind of behavior can be found for the instantaneous free propagator; namely, the 
instantaneous propagator shows the scaling violation even at the tree level (see Appendix B). 
It is the crucial point of the scaling violation that the instantaneous propagator is defined 
as the energy integral of the unequal-time propagator. For a finite temporal lattice spacing, 
the energy integral is limited in the interval [— ri/a T , ir/a T ], and it induces a spurious \p\ 
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dependence on the instantaneous propagator. It leads the discretization errors especially at 
large momenta, as demonstrated in Appendix B, while it disappears only in the Hamiltonian 
limit a T — > 0. Accordingly, the lattice data in the UV region for the instantaneous transverse 
gluon propagator suffer from the discretization errors that cannot be eliminated with the 
cylinder cut or the cone cut we applied. 

One way to circumvent this problem is to exclude the high momentum data from the 
analysis of the propagator. This has been studied in Ref. [16] combined with the matching 
procedure proposed in Ref. [28] in order to find a reasonable value for the available momen- 
tum range. It has been shown that the instantaneous transverse gluon propagator shows 
scaling behavior by restricting the available momentum range. The data on large lattices 
shown in Fig. 2 illustrate that the discretization errors are relatively small in the IR region, 
and this supports the validity of the prescription in Ref. [16]. The another way is to calculate 
the unequal-time propagator _D(j9,p 4 ) and eliminates the p 4 dependence of D(p,p 4 ) as was 
discussed in Ref. [15], and the instantaneous transverse gluon propagator has been shown to 
be multiplicatively renormalizable. 

V. INSTANTANEOUS TEMPORAL GLUON PROPAGATOR ON THE 
ISOTROPIC LATTICE 
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FIG. 3: The dressing function of the instantaneous temporal gluon propagator in physical units at 
j3 = 5.7 and (3 = 6.0. The renormalization point is set to be 2 [GeV] in the left panel and l[GeV] 
in the right panel. 



In this section, we discuss the isotropic lattice result for the temporal gluon propagator. 
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The dressing function of the instantaneous temporal gluon propagator on the isotropic lattice 
at (3=5.7 and 6.0 is drawn in Fig. 3. The dressing function is normalized such that Z 44 (|p| = 
2[GeV]) = 1 in the left panel and Z M (\p\ = l[GeV]) = 1 in the right panel. The dressing 
function Z44 is constant for \p\ at the tree level. In the Gribov-Zwanziger scenario, this is 
expected to diverge in the IR limit resulting in the confining behavior of the color-Coulomb 
potential, which is necessary condition for color confinement in Coulomb gauge QCD. 

Although the scaling violation is clearly visible in Z 44 as in the transverse gluon propa- 
gator, the scaling issue is worse for the temporal gluon propagator than the transverse one. 
In the case of _D tr , the IR data do not much suffer from the discretization errors and they 
almost fall on top of each other by setting the renormalization point to a small momentum. 
Contrastingly, the discrepancy between two data sets of Z44 corresponding to different lat- 
tice spacings remains both in the IR and the UV regions by changing the renormalization 
point (compare the left and the right panels of Fig. 3). Indeed, the scaling violation of .D44 
is not settled by the a-cut method employed in Ref. [16], indicating that the data at small 
momenta may also be affected by the discretization errors. 

Besides the scaling violation, we observe that the dressing function is suppressed in 
the IR region. For the temporal gluon propagator, we expect that its instantaneous part 
behaves as l/|p| 4 in the IR region since it corresponds to the color-Coulomb potential, 
which rises linearly with distance between a quark and an antiquark [4, 5]. As opposed 
to our expectation, the numerical results on the isotropic lattice show that the dressing 
function Z44 decreases with momentum in the IR region for the coarser lattice. 

The unexpected IR suppression of Z44 may stem from an incomplete isolation of the 
instantaneous part in D44. The unequal-time temporal gluon propagator can be decomposed 
into the instantaneous part and the vacuum polarization part as Eq. (10) in the continuum 
theory, and it gives 



On a lattice with finite a T , e is of the order of a T and we would have a 0(a T ) contribution 
from the polarization term in the instantaneous -D44 on the lattice. The unwanted suppres- 
sion at small momenta may originate from such a contribution of the polarization term, 
and furthermore, the energy integral of the polarization term would produce a spurious \p\ 
dependence as for the transverse gluon propagator, which leads the scaling violation. 

Before closing the section, we summarize the points so far: (1) The transverse gluon 




(20) 
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propagator on the isotropic lattice suffers from the discretization errors at large momenta 
and it can be ascribed to the spurious \p\ dependence coming from the temporal lattice 
cutoff in the energy integral defining the instantaneous propagator. (2) The temporal gluon 
propagator is affected both in the IR and the UV regions by the discretization errors, and 
it may stem from the 0(a T ) contribution of the polarization term in the instantaneous 
propagator besides the limited energy integral. 



VI. INSTANTANEOUS TRANSVERSE GLUON 
ANISOTROPIC LATTICES 
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FIG. 4: The instantaneous transverse gluon propagator on the isotropic lattice (top left), on 
the anisotropic lattices with £ = 2 (bottom left), and with £ = 4 (top right). The results for 
the isotropic lattice and the anisotropic lattice with £ = 4 on the large lattice volume are drawn 
together in one figure for direct comparison (bottom right). The propagator is renormalized to 
unity at \p\ = 2 [GeV]. 
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The instantaneous transverse gluon propagators for several anisotropies are drawn in 
Fig. 4. Compared to the isotropic case, the deviations of three curves corresponding to 
different lattice couplings become moderate on the anisotropic lattice with £ = 2. Further 
increase of £ leads to a nice scaling behavior and the data points for £ = 4 almost fall on top 
of one curve. Accordingly, our results on the anisotropic lattices support our expectation 
that scaling violation observed in the instantaneous transverse gluon propagator disappears 
in the limit £ — > oo (a T — > 0). 

In order to investigate how scaling violation is cured by getting close to the Hamiltonian 
limit, we quantify the difference of two curves by the following function, 



x 2 



^ , p f (pQ - gfM y | ( dm - DfHpij | ^ (21) 



=1 x ' i=l x 



Df(pi) represents the measured value of the propagator at the momentum pi and af t i the 
corresponding statistical error. D™ t (p i ) is the estimated value obtained by a cubic spline 
interpolation of the other lattice data set D c . The subscripts f and c label two different data 
sets, meaning "finer" and "coarser" . The summations extend over the data points at which 
the two lattice data sets overlap in the momentum. 

We compute x 2 defined above for the following data sets: 

(1) (24 3 x 48, (3 = 5.80, £ = 2) and (24 3 x 48, (3 = 6.00, £ = 2) 

(2) (24 3 x 96, /3 = 5.75, £ = 4) and (24 3 x 96, (3 = 5.95, £ = 4). 

We note that the physical volumes are very similar; (1) V ~ 4.29 4 [fm 4 ] and V ~ 2.94 4 [fm 4 ], 
and (2) V ~ 4.31 4 [fm 4 ] and V ~ 2.92 4 [fm 4 ], respectively. For each case, we found 

(1) x 2 /iV DF = 867 

(2) rV^DF = 309, 

where N-^p is the number of degrees of freedom of the \ 2 analysis. Decreasing the temporal 
lattice spacing by a factor of 2, xV^df reduces by about a factor 3. We note that the ab- 
solute value of x 2 is unimportant since it depends on the absolute value of the propagator, 
which can take an arbitrarily large (or small) value by multiplicative renormalization. x 2 
defined above make sense only when we compare the x 2 values under the same renormal- 
ization condition and the fixed ratio of the physical volumes of the lattice data sets to be 
analyzed. 
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In the right bottom panel of Fig. 4, the instantaneous transverse gluon propagator on 
the spatial lattice extent L a = 48 is plotted both for the isotropic lattice and the £ = 4 
anisotropic lattice. We observe that the propagator has a maximum at about 500 [MeV] 
irrespective of the lattice coupling and the anisotropy, and it decreases with the momentum 
in the IR region. Accordingly, the turnover of the transverse gluon propagator survives in 
the Hamiltonian limit. 




r[fm] 



FIG. 5: The unrenormalized instantaneous gluon propagator in position space (the correlation 
function of the gauge fields) on £ = 4 anisotropic lattice. The physical volumes are 2.87 4 [fm 4 ] 
(crosses), 2.92 4 [fm 4 ] (left triangles), 4.31 [fm 4 ] (down triangles), 5.74 [fm ] (up triangles), and 
8.61 4 [fm 4 ] (circles). 

The confined behavior of the gluon propagator can be seen more directly in position space. 
The unrenormalized correlation function of the transverse gauge fields on the anisotropic 
lattice £ = 4 is depicted in Fig. 5. On small lattices (crosses and left triangles in the figure), 
the correlation function is positive and a decreasing function of the distance. As the physical 
volume increases, D tv (r) develops a dip at about r = 1 [fm] and becomes negative around 
the dip. On the largest volume, the correlation function quickly decreases with distance at 
small r and becomes negative in the range r ~ 1 — 2 [fm] and then vanishes at large distances. 
It means that the gauge fields have no correlation over the hadronic scale in sharp contrast 
to the massless particles for which the correlation function behaves as 1/r 2 . We notice that 
the results at the approximately fixed physical volume nicely agree with each other; namely, 
lattice discretization errors do not seriously affect the correlation function on the anisotropic 
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lattice. 



VII. POWER LAW FITTING OF D tr AT SMALL AND LARGE MOMENTA 

We note that the momentum dependence of the instantaneous transverse gluon prop- 
agator on the anisotropic lattices in the UV region significantly differs from that on the 
isotropic lattice (see the bottom right panel in Fig. 4). On the isotropic lattice, the slope 
at high momenta gets smaller as the lattice couplings increase. On the anisotropic lattices, 
the slope decreases further. In order to investigate the UV behavior of the transverse gluon 
propagator, we make the power law ansatz, 

D "® = h*3" (22) 

and fit this ansatz to the data on 24 3 x L T lattices in the momentum range \p\ > 6 [GeV]. The 
fitted parameters are given in Table II. We see that 7^} v gets small values as the anisotropy 
increases, and that for £ = 4 takes about one third of that on the isotropic lattice. 

As we have discussed in Sec. IV, the inverse of the instantaneous propagator can be 
interpreted as the energy dispersion relation. Since gluons are expected to behave as free 
massless particles at sufficiently large momentum due to asymptotic freedom, this interpre- 
tation can be allowed if the instantaneous transverse gluon propagator behaves as l/|p| in 
the UV region, which corresponds to the null UV exponent. Therefore, the decrease of 7^} v 
with increasing £ is consistent with the asymptotic free field behavior of the gluon fields in 
the continuum limit 2 . 



At least, the linear extrapolation of 7^ to the Hamiltonian limit gives 7gj (£ — > 00) = 0.0852(201) 
and x 2 /-/Vdf = 11.1, that is, the value of 7^ v decreases but remains Unite in this analysis with three 
anisotropies, although the linear extrapolation may be too naive and inadequate, and a careful study of 
taking the Hamiltonian limit is needed. 
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TABLE II: The result of the UV power law fitting of D tr in the momentum range \p\ > 6 [GeV]. 

(L a ,L T ,£,P) <h 7^ x 2 /iVDF 

(24, 24, 1, 6.0) 4.89(14) 0.920(15) 1.34 

(24, 48, 2, 6.1) 3.83(25) 0.617(35) 0.86 

(24, 96, 4, 6.1) 2.49(7) 0.282(15) 0.68 



In order to explore the IR behavior of the instantaneous transverse gluon propagator, we 
make the power law ansatz, 

D tr (p) = d 2 \p\^, (23) 

in the IR region. The fitted parameters are listed in Table III. 

Although the fitting becomes worse and the IR exponent 7*f becomes small as the max- 
imum momentum of the fitting range increases, 7*1 takes a positive value in all the cases. 
For the anisotropic lattice, we need much larger lattices to extract the IR exponent with 
an acceptable \ 2 value. In both the isotropic and the anisotropic cases, our result of the 
IR fitting predicts the vanishing transverse gluon propagator at zero momentum. Given the 
fact that the fitted values of 7^ increases with decreasing the maximum momentum of the 
fitting range, the 7*f values listed in Table III would give a lower bound for the IR exponent 
of the transverse gluon propagator in the Coulomb gauge. 

We note that fitting the transverse gluon propagator with the Gribov-type ansatz, Eq. (2), 
which successfully reproduced the lattice data in Ref. [15], did not work in our lattice data 
for both the isotropic and the anisotropic cases. At the same time, the peak position of 
the transverse gluon propagator differs between our results and the results in Ref. [15] ; it is 
about 500 [MeV] in our case and the fitting analysis in Ref. [15] gives 880(10) [MeV], which 
is rather close to the peak position of the dressing function Z tr in our results (see the right 
panel of Fig. 2). It deserves further study to clarify whether this discrepancy comes from 
the gauge group [SU(3) in this study and SU(2) in Ref. [15]] or the adopted prescriptions to 
circumvent scaling violation. 
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TABLE III: The result of the IR power law fitting of D tr . \p ma x\ represents the maximum 
momentum of the fitting range. 

0.27 6.64(17) 0.311(17) 0.274 
( 48 ~ 64, 1, 5.7 ) 0.30 6.24(9) 0.271(10) 4.48 
0.33 5.98(6) 0.241(8) 7.82 
( 48, 4, 5.75 ) 0.45 4.09(2) 0.174(4) 94.3 



VIII. INSTANTANEOUS TEMPORAL GLUON PROPAGATOR ON 
ANISOTROPIC LATTICES 

We next discuss the instantaneous temporal gluon propagator, which is related to the 
color- Coulomb potential. The dressing function of the time-time component of the gluon 
propagator is shown in Fig. 6 for the isotropic lattice and the anisotropic lattices. On the 
anisotropic lattices, the dressing function shows much better scaling behavior than that on 
the isotropic lattice. Although the small deviation can be seen both in the IR and UV 
region, one can expect that the scaling behavior is completely recovered in the Hamiltonian 
limit. 

Moreover, we find that the IR behavior of Z44 on the anisotropic lattice is different from 
that on the isotropic lattice (see bottom right panel). For the isotropic case, we see that the 
dressing function increases slowly with decreasing the momentum, and bends down on coarse 
lattice. By contrast, Z 44 continues to rise on the anisotropic lattice even for the coarsest 
lattice data (/3 = 5.75), and Z44 at available smallest momentum for the anisotropic case 
is about 10 times larger than that for the isotropic case. We note that the spatial lattice 
spacing for (£,/?) = (4,5.75) is larger than that for (£,/?) = (1,5.70). This implies that 
Z44 is very sensitive to the discretization effects, and taking the Hamiltonian limit is crucial 
to cure scaling violation for the temporal gluon propagator and to explore the genuine IR 
divergent behavior in the Coulomb gauge QCD. 
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FIG. 6: The dressing function of the instantaneous temporal gluon propagator on the isotropic 
lattice (top left), on the anisotropic lattice with £ = 2 (bottom left), and that with £ = 4 (top 
right). The results for the isotropic lattice and the anisotropic lattice with £ = 4 on the large 
lattice volume are drawn together in one figure for direct comparison (bottom right). The dressing 
function is renormalized to unity at |p| = 2 [GeV]. 

IX. POWER LAW FITTING OF Z 44 AT SMALL AND LARGE MOMENTA 



The asymptotic form of the color-Coulomb potential, the instantaneous part of the tem- 
poral gluon propagator, is given by 



|piv 44 (bi; 



%0 



(24) 



26 ln(|p1/A C oui)' 

where xo = 12/11 and b = 1 1 /16tt 2 for SU(3) pure Yang-Mills theory, and Ac ou i is a finite 
QCD mass scale [29]. We here fit the data for Z 44 with 



Zu(\p\) = Z 



x 



26 ln(|p1/A C oui) 



(25) 
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in the momentum range \p\ > 6 [GeV]. The fitted parameters Z and Acoui are given in 
Table IV with x 2 I^df- We find that the fitted parameter A Cou i takes an unacceptably 
small value as a QCD mass scale, although the x 2 value is reasonable for the isotropic 
lattice. Acoui increases with the anisotropy and the fitted value Ac ou i = 0.8345(85) [GeV] 
for £ = 4 is the order of Aqcd- However, the numerical data for the Aqcd parameter are still 
unstable under the increase of the anisotropy and Aqcd could change by further varying £. 
It can be stated that the anisotropy should be greater than 4 otherwise we have physically 
inadequate results, although it is difficult to estimate the Hamiltonian limit of Aq CD . 

TABLE IV: The result of the UV logarithmic law fitting of Z44 in the momentum range \p\ > 6 
[GeV]. 

(Z4xL r ,g,/3) Z A C oui[GeV] X 2 /iV DF 
(24 4 , 1, 6.0) 0.648(97) 0.000968(1282) 0.406 
(24 3 x 48, 2, 6.1) 0.310(80) 0.0385(507) 1.20 
(24 3 x 96, 4, 6.1) 0.0842(41) 0.845(85) 2.22 



In the IR region, the color-Coulomb potential is expected to behave as l/|p| 2 , which 
gives a linearly rising potential in position space. Indeed, the lattice simulations of the 
color-Coulomb potential obtained from the correlator of the partial Polyakov line revealed 
such a linearity of the color-Coulomb potential at large distances [4, 5]. Therefore, we fit 
the dressing function Z 44 with the power law form, 

Z «dPl) = < 26 > 

Assuming that the contribution of the vacuum polarization term to the instantaneous tem- 
poral gluon propagator is negligible, we expect 744 = 2 giving a linearly rising color-Coulomb 
potential. 

The fitted results are given in Table V. We observe that x 2 is extraordinary large and the 
fitting does not work. Even though the scaling violation becomes moderate as the anisotropy 
increases, we still have discretization effects which are visible in Fig. 6; the curves associated 
with different lattice couplings deviate from each other in the IR region and the slope gets 
steeper as the lattice spacing decreases. Therefore, we are still not close to the Hamiltonian 
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limit where the scaling behavior is observed, and we cannot extract the color-Coulomb string 
tension from ^44 ( |p| ) and compare it with that obtained from the link-link correlator. 

TABLE V: The result of the IR power law fitting of Z44. |p max | represents the maximum momen- 
tum of the fitting range. 
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X. INSTANTANEOUS D 44 IN POSITION SPACE AND THE COLOR-COULOMB 
POTENTIAL 




FIG. 7: The instantaneous temporal gluon propagator in position space on the isotropic lattice 
(left) and the £ = 4 anisotropic lattice (right). 

Since the instantaneous part of the temporal gluon propagator corresponds to the color- 
Coulomb potential, 

D M (x, t) = V c (x)S(t) + P(x, t), (27) 

the color-Coulomb potential in the color-singlet channel is given by measuring the instanta- 
neous -D44 in position space, 

Vr^(x) = -^a T D M (x). (28) 
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Here the factor —4/3 is the Casimir invariant in the fundamental representation in the SU(3) 
gauge group. The instantaneous temporal gluon propagator in position space is drawn in 
Fig. 7 for the isotropic and the £ = 4 anisotropic lattices. We observe that the behavior 
of .D44 notably changes by increasing the anisotropy. decreases almost linearly with 
distance in the range 3 < r < 7 [fm] for the anisotropic lattice while it does not on the 
isotropic lattice. The linear decrease of -D44 means that V c is a linearly rising potential with 
distance, which is consistent with the lattice calculations of V c from the correlator of the 
partial Polyakov line. On the isotropic lattice, there are unavoidable contributions from 
the polarization term in the instantaneous propagator on the lattice, and it is indispensable 
to carry out the lattice simulations with small temporal spacing to extract V c from the 
instantaneous -D44. 

We fit the data with the function 

V c sinslct (r) = C + a c r, (29) 

in the range 0.4 < r < 0.7 [fm] and extract the color-Coulomb string tension a c . The 
fitted result is given in Table VI. We see that the fitting is extremely worse for the isotropic 
lattice, and xV^df approaches 0(1) on the anisotropic lattice. Besides a large \ 2 value, 
the color-Coulomb string tension on the isotropic lattice is smaller than the Wilson string 
tension violating the Zwanziger's inequality [8]. The extracted a c increases with decreasing 
the temporal lattice spacing and satisfies the Zwanziger's inequality on the finer lattices. 
Although a c from the instantaneous D44 is still smaller than that from the correlator of 
the partial Polyakov loop, there is little doubt that a c does not saturate the Zwanziger's 
inequality and a c is larger than the string tension of the Wilson static potential. 

It may be surprising that we can extract a c with reasonable y 2 values on finer anisotropic 
lattices. We have seen in the last section that the power law fitting of Z44 does not work, 
even though the extracted value of 744 is close to 2, which gives a linearly rising color- 
Coulomb potential. On the anisotropic lattice, the largest lattice volume is 4.67 4 [fm 4 ] 
(circles in the right panel of Fig. 7) and we naively expect that the finite volume effects on 
the instantaneous -D44 are not so serious in the range r < 2 [fm]. However, Fig. 7 shows that 
-D44 does not show a linearly decreasing behavior at r > 1 [fm]. This indicates that we still 
need to approach the Hamiltonian limit to extract the color-Coulomb potential from the 
instantaneous temporal gluon propagator and to compare the color-Coulomb string tension 
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with that from the correlator of the partial Polyakov loop. 

TABLE VI: The result of the fitting of V c in the range 0.4 < r < 0.7 [fin]. 



(L a ,L T ,£,(3) C [GeV] [MeV] X 2 /iV DF 
(48, 48, I, 6.00) -0.4971(4) 289.9(3) 39.1 

(48, 192, 4, 5.75) -1.262(2) 423.2(9) 11.3 

(48, 192, 4, 5.95) -1.995(5) 515.7(18) 1.26 

(48, 192, 4, 6.10) -2.634(5) 580.0(16) 0.828 



XI. CONCLUSIONS 

It is a central interest to explore the IR behavior of the gluon propagator to reveal the 
confinement mechanism in QCD. In the Coulomb gauge, the instantaneous gluon propagator 
suffers from significant discretization errors on isotropic lattices. In this paper, we calcu- 
lated the transverse and the temporal components of the instantaneous gluon propagator 
on isotropic and anisotropic lattices and studied the scaling behavior. 

We find that the transverse gluon propagator shows a nice scaling behavior on the 
anisotropic lattices, and scaling violation observed on the isotropic lattice almost disap- 
pear on £ = 4 anisotropic lattice. It is natural to expect that the perfect scaling behavior 
can be seen and the multiplicative renormalizability holds in the Hamiltonian limit a T — > 0. 
In the IR region, the transverse gluon propagator is strongly suppressed on both the isotropic 
and the anisotropic lattices and shows the turnover at about 500 [MeV] in both cases. 

We also calculated the transverse gluon propagator in position space, and it shows that 
the correlation function quickly decreases at small distances and becomes negative in the 
range r ~ 1 — 2 [fm], and vanishes at large distances. This means that the gluon fields have 
no correlation beyond the hadronic scale, and it is consistent with the fact that the gluons 
are confined in the hadrons (glueballs). 

The power law fitting of the transverse gluon propagator exhibits that the UV exponent 
7^{ v decreases with increasing the anisotropy. This supports the expectation that the gluons 
behave as free massless particles at sufficiently large momentum in the continuum due to 
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asymptotic freedom. 

In order to explore the IR behavior of the transverse gluon propagator, we fitted the 
data with the power law and found that the extracted value of the IR exponent is positive 
on the isotropic and anisotropic lattices. However, our lattice is still small to extract a 
reliable value of the IR exponent on the anisotropic lattices. Fitting the transverse gluon 
propagator with the Gribov-type ansatz did not work in our lattice data, which successfully 
describe the SU(2) lattice data employing different method to cure scaling violation [15], 
and it deserves further study whether this is due to the difference of the gauge group or the 
adopted prescriptions to remedy scaling violation. 

For the temporal gluon propagator, scaling violation was observed both in the IR and 
UV regions on the isotropic lattice in contrast to the transverse propagator, for which lattice 
data do not suffer from discretization errors in the IR region. Our results show that the 
lattice data on the anisotropic lattices show much better scaling behavior than that on 
the isotropic lattice, although the discretization errors are still seen in the IR and the UV 
regions. 

We observed that the time-time gluon propagator on the anisotropic lattice is much 
more enhanced in the IR region compared to that on the isotropic lattice. The turnover 
observed on the coarsest lattice data on isotropic lattice disappears on the finer lattices, and 
Z44 monotonically increases with decreasing the momentum on £ = 4 anisotropic lattice. 
However, the IR power law fitting did not work for the dressing function, although the fitted 
value 744=1.84(1) on the finest lattice is close to the expected value for the linearly rising 
behavior of the color-Coulomb potential. We again need larger lattices as the transverse 
gluon propagator. 

The logarithmic law fitting of the dressing function of the temporal propagator in the 
UV region indicates that the extracted A Cou i is unacceptably small on the isotropic lattice 
and it takes 0.845(85) [GeV] on the finest anisotropic lattice, which is the order of A QCD . 

In position space, the linearly decreasing behavior can be seen on the anisotropic lattice, 
and the color-Coulomb string tension obeys the Zwanziger's inequality on finer anisotropic 
lattices. However, the extracted value is still smaller than that obtained from the correlator 
of the partial Polyakov line. Thus, we can say that although the scaling violation is softened 
by decreasing the temporal lattice spacing, the instantaneous temporal gluon propagator 
receives a contribution from the polarization term and it is difficult to extract the color- 
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Coulomb potential from the instantaneous -D44 for finite a T . 
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Appendix A: Tolerance of the gauge fixing and the IR suppression of D tr 

The tolerance of the gauge fixing is crucial for the IR suppression of the instantaneous 
transverse gluon propagator. In the iterative gauge fixing on the lattice, we stop the gauge 
transformation if the violation of the transversality becomes less than some small number e; 

d = {N ^ l)L ^ ^?{xM<e. (Al) 

We set e = 10 -14 in our calculations. In this appendix, we examine the behavior of the 
instantaneous transverse gluon propagator by varying e on the isotropic lattice. 




12 3 4 

pa 



FIG. 8: The instantaneous transverse gluon propagator in lattice units on 56 4 lattice at [3 = 5.7 
for various e. The cylinder cut is applied, but the cone cut is not. 

Figure 8 shows D tr (p) in lattice units on 56 4 lattice at /3 = 5.7 for various e. For all the 
cases, the measurements are done for 20 configurations and the cylinder cut is applied. We 
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observe that the propagator does not exhibit a turnover for e = 10~ 2 and the curve deviates 
from that for e = 10~ 14 in the region pa < 1.5. Decreasing e by the factor 100 diminishes 
the propagator at small momenta and D tr coincides within the statistical errors in the range 
1 < pa. A further decrease of e suppresses D tT in the IR region and the transverse propagator 
shows a clear turnover. Although the data for e = 10~ 8 are not depicted in Fig. 8, they fall 
on top of the data for e = 10 -14 within the error bars. Therefore, e should be sufficiently 
small to explore the IR behavior of the gluon propagator and 10 -14 is small enough. 

Appendix B: Scaling violation in the free field case 

In this appendix, we discuss scaling violation of the instantaneous propagator in the free 
field case. A more general case was discussed in Ref. [15]. 

The instantaneous transverse propagator in the continuum is given by integrating the 
unequal-time propagator over p 4 , 



In the free field case (i.e., at the zeroth order in the coupling), the four-dimensional propa- 
gator is Df iee (p, p^ = l/(|p| 2 + pi) and the instantaneous propagator is 



On a lattice with a finite temporal spacing a T , the p^ integral is limited within the range 
—n/a T < p4 < Ti /a T . Thus, the instantaneous propagator is given by 




(Bl) 




1 



(B2) 



2bT 




(B3) 



For a finite temporal lattice spacing, we have an extra factor 




(B4) 
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where p = \p\a a is the spatial momentum in lattice units, and £ is the anisotropy, £ = a a /a T . 
This extra factor approaches unity in the Hamiltonian limit a T — > 0. 




i 

p [GeV] 



FIG. 9: Df ree (\p\, a T ) is plotted as a function of p for various temporal lattice spacings. 

The free instantaneous propagator for various temporal lattice spacings is illustrated in 
Fig. 9. We observe that the propagator starts to deviate from that in the Hamiltonian limit 
at larger momenta as a T decreases. Accordingly, scaling violation is observed even in the 
free field case and it goes away in the Hamiltonian limit. If we impose some renormalization 
condition, Df vee (\p\ = l[GeV]) = 1 [GeV -1 ] for instance, the curves in Fig. 9 cross at the 
renormalization point and deviate from each other in both the small and the large momentum 
regions, which is the situation we encountered in the lattice simulations. From this simple 
exercise, we expect that scaling violation for the instantaneous gluon propagator would be 
moderate for small temporal lattice spacing. 
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